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Abstract 



We study the spin-boson model with a sub-Ohmic bath using infinitesimal unitary transforma- 
tions. Contrary to some results reported in the literature we find a zero temperature transition 
from an untrapped state for smah coupling to a trapped state for strong coupling. We obtain an 
exphcit expression for the renormalized level spacing as a function of the bare papameters of the 
system. Furthermore we show that typical dynamical equihbrium correlation functions exhibit 
an algebaric decay at zero temperature. 
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The problem of a small system coupled to environmental degrees of freedom plays an impor- 
tant role in various fields of physics. The spin-boson model is one important example in this 
class of models. Here a two level system, i.e. a pseudo spin, is coupled to environmental degrees 
of freedom modeled by non-interacting harmonic oscillators, i.e. bosons. It is defined by the 
Hamiltonian 
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H = --Aoa^ + -a,Y,Xk{bl + bk) + J2iOkbibk + Eo . (1) 
^ ^ k 

Aq is the level splitting of the two level system. The couplings and and the bath frequencies 
LOk enter only via the combination J{uj) = ^k'^'i^i^ ~ ^k)- This model is a paradigm in the 
theory of dissipative quantum systems. It has a wide range of applications. For instance it has 
been used to study dissipative tunneling. In this case az plays the role of the coordinate of the 
tunneling particle, Aq is the bare tunneling frequency. We will use this language in this letter. 
The model has been investigated by various methods. For a review of applications, methods and 
results we refer to 

One usually assumes that the spectral function J{uj) behaves as uj^ with s > at low frequen- 
cies. Furthermore J(w) should contain some high frequency cutoff lOc- The case s = 1 is called 
the Ohmic bath, it has been widely studied and is well understood. As a function of the coupling 
strength the model shows a zero temperature transition from a trapped to an untrapped state. 
This has been shown rigorously by Spohn and Diimcke Q. They mapped the spin-boson model 
to an Ising model on with a long-range interaction by integrating out the bosonic degrees of 
freedom. The corresponding Ising model shows a phase transition for < s < 1. The point s = 1 
is the phase boundary, for s > 1 the Ising model never shows long-range order. The ordered 
phase of the Ising model corresponds to the trapped state in the spin-boson model. Thus in 
the super-Ohmic case s > 1 the particle is never localized, whereas in the Ohmic case one has 
a transition from an untrapped state for small coupling to a trapped state for sufficiently large 
coupling. 

In the literature about the spin-boson problem exists some confusion about the existence 
of such a transition for the sub-Ohmic bath s < 1 as within some approximation schemes the 
particle is localized for any nonzero sub-Ohmic coupling at zero temperature 0. In contrast 
the corresponding results for the Ising model on ?R. indicate a transition as a function of the 
coupling strength Q. However, for s < 1 these results for the Ising model have not been carried 
over to the spin-boson model rigorously. Still, based on the findings of Spohn and Diimcke 
it would appear surprising not to find a corresponding transition in the spin-boson model. It 
should be noted that the sub-Ohmic case is of less physical interest as compared to Ohmic or 
super-Ohmic coupling, still it contributes to the general understanding of the parameter space 
of the spin-boson problem. Also one can learn something about possible problems of adiabatic 
renormalization schemes as we will see. 

In Ref. 0] we recently studied the spin-boson model with an Ohmic bath using a new method 
proposed by Wegner This method simplifies a given Hamiltonian using infinitesimal unitary 
transformations. The advantage of infinitesimal unitary transformations is that one is able to 
eliminate couplings between states with a large energy difference first, couplings between states 
with smaller energy differences later. In this manner infinitesimal unitary transformations auto- 
matically seperate different energy scales. A continuous unitary transformation can be written 
as a differential equation for the Hamiltonian or equivalently as a set of coupled differential 
equations for the coupling constants. In the case of a dissipative quantum system, the goal of 
the transformation is to decouple the system from the bath. This yields an effective Hamiltonian 
for the system, which can be used to calculate e.g. dynamical quantities. In the case of the 
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spin-boson model the effective Hamiltonian is only a two-state system that can be described 
by an energy shift and a renormalized level spacing. The advantage of this method is that the 
unitarity of the time evolution of system plus bath is not destroyed. One is able to obtain a 
systematic expansion for the long-time dynamics and the static behaviour of the system. 

In the present letter we also study the sub~Ohmic case in some detail using this method. 
Our results will turn out to be in agreement with the rigorous results in Ref. Q. In Ref. 
we derived an equation for the renormalized tunneling frequency Aqo as a function of the bare 
tunneling frequency Aq and the spectral function J{uj). We sketch the derivation of this equation 
and use it to calculate the renormalized tunneling frequency in the sub-Ohmic case. Whereas 
Aoo is a continuous function of the coupling strength for the Ohmic bath, the transition seems 
to be discontinuous in the sub-Ohmic case. We calculate an explicit expression for the critical 
value of the bare tunneling frequency for which the transition occurs. Furthermore we apply 
some general results from Ref. concerning time-dependent equilibrium correlation functions 
to the sub-Ohmic case. 

A continuous unitary transformation of the Hamiltonian can be written as a differential 
equation for the Hamiltonian 

^^ = [^{l),H{l)] (2) 

with the initial Hamiltonian H{Q) = H given in Eq. (|l|). r] is the generator of the transformation 
and depends on the flow parameter i. r] is an antihermitian operator and we choose 

i 1 

^ = 2^^ ^ %,3/(&fc + bi) - -(Tz mAbk - bl) . (3) 

k k 

r]]. y and rjk^z are ^-dependent parameters. Performing the transformation the Hamiltonian be- 
comes ^-dependent as well. In addition to the coupling of cj^ to the bosonic modes, new couplings 
are generated. One of these new terms is of the form ay Yl,k Ck{bk — b]^). It can be avoided by a 
suitable choice of the parameters ij^^y and ij^^z- The other terms describe a coupling of the spin 
to more than one bosonic mode. In Ref. Q| we showed that within a very good approximation 
one can restict oneself to a linear coupling of the spin to the bath. Thus one can neglect the 
normal ordered part of the higher interactions. As a consequence the Hamiltonian has the form 
given in (|l|) for any i, but now the parameters and A depend on £. Using 

Vk,y = -^>^kf{^kJ), Vk,z = ^k^kfi^kJ) (4) 

the £ dependence of the parameters is given by the flow equations 

dXk 
di 



-Xk{ujl-A^)f{uk,e), (5) 



^ = -Aj2>^lfi^kJ){2nk + l) , (6) 

k 

^ = -\Y-^kXlfiu^k,i), (7) 

k 

where =< b^bk >= {exp{(3ujk) — 1)^^- For details of the calculation we refer to Ref. |^. 
Introducing the ^-dependent spectral function J{lo,£) = Ylik ^'k(^)^(.^ — ^k) the flow equations 
can be written as 

^ = - A / dujJiio, i)fiu;, i) (2n(^) + 1) , (8) 
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^^^ = -2{u;'-A')f{u;,e)J{u;,i), (9) 

with the occupation number n{uj) = [e^^ — 1)"^. These two equations can be combined into a 
single equation 

dlnA _ 1 r dJ{uj,l)2n{uj) + l 

o / «^ -71 Xq- ■ V^^i 



d£ 2 J de cj2-A2 

Eq. ( |lO|) can be solved approximately by replacing A{£) by Aqo on the right hand side as explained 
in Refs. 0, This yields the self-consistency condition 

\n— = -- duJ{uj,0)— -2-. 11 

Ao 2 J u;^ - A^ 

The integral on the right hand side contains a singularity, it is to be interpreted as a principal 
value integral. This equation is applicable for any spectral function J(cj, 0). The only restriction 
is that the coupling must not be too strong, since then processes involving more than one bosonic 
mode become important. It can be shown that if such proccesses are included initially in the 
Hamiltonian, an equation of the form (|Tl|) is obtained as well, but with a modified spectral 
function J(ti;,0). Thus the main problem is not the occurence of multi-boson processes, but the 
correlations between the bosonic modes which are induced if the coupling is too strong. 

A general feature of equation ( [TTI) is that slow modes lo < A^o tend to increase the renormal- 
ized level splitting whereas fast modes tend to decrease it. This effect can already be observed 
in a perturbational treatment. The balance of these effects can be essential when one treats a 
problem with a discrete eigenvalue embedded in a contiuum (compare for example Ref. for 
the case of the Anderson impurity model). For the spin-boson problem this balance becomes 
particularly important for the sub-Ohmic bath as will be discussed in the sequel. 

In Ref. Q the function f{uj,i) in Eq. (^) was chosen as 

/(.,f) = ^. (12) 

With this choice the coupling function J{uj,£) decays monotously as a function of i for all values 
of CO as depicted in Fig. 1. For the Ohmic bath J{cu,0) = 2aLo9{uJc — ^) the numerical solution 
of the differential equations with this choice of f{Lo,i) turned out to be in good agreement with 



the self-consistency equation (pTf). 



For the non-Ohmic bath the coupling is parametrized as 

J{uj,0) = K^''uj'e{uJc-uj) (13) 

with a coupling constant K with the dimension energy. The explicit form of the cutoff function 
is not important. In contrast to the Ohmic bath, the non-Ohmic bath contains a second energy 
scale given by the coupling constant K. Whereas lOc is usually assumed to be much larger than 
Aq, this new energy scale can be of the order of Aq. For the super-Ohmic case s > 1 and T = 0, 
we can neglect A in the integral on the right hand side of (^) or A^o in the integral on the 



right hand side of (11) for Aq <C uJc since the integral is not infrared divergent. This yields the 
well-known result 

A~ = Aoexp(-^^). (14) 

Comparison with a numerical solution of the flow equations shows that this is a very good 
approximation to the solution of the differential equations for Aq ^ lOc- For Aq > uJc the flow 
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equations yield a solution A(^) that increases with £ and Aqo is larger than cOc- In this limit 
the formula ( p!^ ) is not valid. For the super-Ohmic bath the particle is never localized, there is 
no transition as in the case of Ohmic dissipation. Our value of Aqo in (|T^ ) is consistent with 
adiabatic renormalization [0. 

We are interested in the sub-Ohmic case. For s < 1 Aqo can be calculated in the limit 
LOc ^ oo since the integral in (|ll|) is not ultraviolet divergent. Then (11) yields 




Ao = Aoo exp ( ^ ) cot(-(l-s))) (15) 
If the condition 



def ,^/4tan(f(l-s))^ 



1 

s-l 



Ao > A, = i^ r—y^ (16) 

V e^(l-s) J 
is satisfied, the solution for A^o is positive and satisfies 



'4tan(f(l-s))^ 



1 

s-l 



Otherwise Aqo = 0. For the sub-Ohmic bath we obtain a transition from an untrapped state 
for sufficiently large Aq as compared to K to a trapped state for small Aq. In contrast to the 
transition in the Ohmic case, which depends only on the dimensionless coupling strength a, the 
present transition depends on the ratio between Aq and K. It occurs for Aq = Ac. Let us remark 
that Ac oc is a simple consequence of the fact that K can be used to fix the energy scale. 
This has also been obtained by Spohn and Diimcke using a variational calculation [^. From 
the self-consistency condition (|Tl|) we therefore conclude that the transition is discontinuous 
since we obtain Aqo = Acexp(^^) > for Aq = Ac. As a cautionary remark we would like to 
emphasize that the nature of the transition is difficult to investigate from the numerical solution 
of the original differential equation. Discontinuous behaviour is only possible for the asymptotic 
solution and hence extremely sensitive to accumulation of errors. But away from the transition 
point we find good agreement with the self-consistency condition and discontinuous behaviour 
at the transition point is more compatible with the numerics. In any case the existence of a 
transition from an untrapped to a trapped state for increasing coupling K is certain. 

Using (|l^) it is a straightforward problem to study Aqo or the critical value Ac as a function 
of ujc- The results (|l5| ) and (|l^ ) are only valid if iOc is large compared to all other energy scales. 
For smaller values of uJc our result Ac (x K no longer holds. 

It is also possible to treat the problem for finite temperature T > 0. The reasoning is identical 
to Ref. Q. For finite temperature Eq. ( p4| ) goes over into 

Ao = A^ exp (^) cot(^(l - .))) F (^, |) (18) 
where 



For small values of T < A^o one finds F < 1. Therefore Aao{T) increases with T for small T. 
For large T, F becomes large and Aao{T) decreases and finally tends to zero for some critical 
temperatue Tc- This behaviour is similar to the Ohmic case |Q. The transition from a finite value 
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of Aoo at low temperatures to Aqo = at high temperatures can be interpreted as a transition 
from coherent to completely incoherent motion of the tunneling particle. 

Equilibrium correlation functions at T = can also be calculated as in Ref. Q. We consider 
the autocorrelation function 

1 1 

C{t) = ^{(^z{t)aM + a,(0)cj,(t))^ = -Tr(pegk^(t)a,(0) + a,(0)a,(t)]) , (20) 
where 

Peg = exp{-PH)/Tr(^expi-(3H)) . (21) 

Since the Hamiltonian has a simple form for ^ = oo, we want to calculate the expectation values 
in this limit. The advantage is that for £ = oo the equation of motion for operators becomes very 
simple and can be solved explicitly. Therefore we do not destroy the unitarity of time evolution 
in contrast to many other approximations. Furthermore it is very easy to perform the thermal 
average. The price we have to pay is that we have to calculate the operator of interest at £ = oo. 
To do this we must solve the flow equation 

^^ = [rj{£),a.{£)] , a,(0) = a. . (22) 

Clearly the expression for (Tz{£) is rather complicated. Fortunately it turns out that one can 
obtain the correct long-time dynamics with the ansatz 

az{£) = h{£)az + a^Y. Xkimk + hi) , (23) 

k 

i.e. terms containing more than one bosonic operator are neglected. We analyzed the flow 
equation for h{£) and Xfc(-^) in Q- I^i the following we report the general results obtained in 
and apply these to the sub-Ohmic bath. The Fourier-transform of C{t) = ^ C{u!) e"*"^* du 
can be written in the form 

C{uj) = l-^d{uj-Aoo) + - ^ hAf{uj,£)^Jiu;,£)d£j e{oj) + (co ^ -co) 




^5(a;-Aoo) + vr / ^ ^/ ' d£\ e{uj) + {u: ^ -u:) 




A2 d£ 

(24) 

If Aoo = one finds /i(oo) 7^ as expected. Hence the particle is localized and C{t) tends to 



a constant for large t. If Aqo > one has h{oo) = and (24) is correct to leading order in a;, 
i.e. corrections are at least of order J(a;,0)^ for small w. For uj <C Aqo, J{uj,£) tends to zero 
rapidly. Therefore we can replace h and A in the integral by their initial values. This shows that 
C{uj) oc J{uj) oc uj'^ for small uj and consequently 

C{t) oc t-'-^ (25) 

for large times. This long-time behaviour of the equilibrium correlation function holds for gen- 
eral s if the particle is not localized. Therefore the asymptotic dynamics is a power-law relaxation 
in the disordered phase at zero temperature. Let us mention that Monte-Carlo simulations for 
this correlation function in the sub-Ohmic case reported in |Q] indicate an algebraic decay. 

As already mentioned in the beginning of this letter, there are various claims in the literature 
that the particle is always localized in the sub-Ohmic case for zero temperature |l|, |2| in con- 
tradiction to our results reported here. It is interesting to analyze the reasoning leading to this 
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conclusion. One important method for treating the spin-boson problem is the "non-interacting 
blip approximation" (NIBA) In order to apply this scheme, the typical blip width must be 
much smaller than the typical distance between blips. For the sub-Ohmic bath this is equivalent 
to Aq <C -ftT [||]. Then the NIBA predicts localization, which is obviously in agreement with 
our condition (16). When Aq and K are of the same order of magnitude, the NIBA cannot be 



justified irrespective of the value of lOc- Anyway the NIBA is unreliable for studying the long- 
time behaviour of equilibrium correlation functions as it gives an exponential decay in Eq. ( p5| ) 
instead of an algebraic decay |l|] . In this sense we agree with the conclusions of Ref . |^] . For large 
times and T = our method of continuous unitary transformations is a simpler and systematic 
approximation scheme. 

Another tool for analyzing the spin-boson model is adiabatic renormalization. This relies 
on successively integrating out high frequencies as compared to the tunneling frequency A. If 
one integrates out the frequencies between some lower cutoff uji » Aq and Wc one obtains a 
renormalized tunneling frequency A'{uji) in the zeroth-order adiabatic approximation (for details 
of the calculation see Ref. [Q]) 

A'{u;i) = Aoexp(-l r ^ duj\ (26) 



. 2 7^, w2 

Choosing uji = p Aq with some number p large compared to unity, the above procedure can be 
iterated: Since A'{ijJi) < Aq one next integrates out the frequencies between uj'i = pA'{uJi) and 
uji and so on. This procedure is depicted in Fig. 2. The iteration converges to a finite value of 
the tunneling frequency if a nonzero solution of 

Aoo = Ao exp (-^[^^ A-i^ (27) 

exists (we have taken the limit ujc oo). But this depends on the ratio K/p. Above a certain 
critical value only the trivial solution A^q = exists. Strictly speaking this renders the adiabatic 
renormalization scheme useless in the sub-Ohmic case as localization or non-localization cannot 
depend on the arbitrary number p. On the other hand one easily checks that p plays no similar 
role in the Ohmic or super-Ohmic case. 

In contrast in the method employed in this letter the modes are decoupled in a continuous 
manner as depicted in Fig. 1. Here corresponds to the energy difference that is just being 

decoupled for a certain value of i. The fundamentally important separation of energy scales in 
the adiabatic renormalization scheme (or more generally in the "poor man's" scaling approach) 
is with respect to the total energy of the modes as measured from some fixed zero point of energy, 
whereas our procedure separates any two modes with respect to their specific energy difference. 
To deal with large energy differences first and then proceeding to smaller energy differences 
is intuitively the safe order of things, e.g. this is also one of the basic concepts underlying 
Wilson's numerical renormalization approach |^]. From this viewpoint one can also understand 
why adiabatic renormalization works for the Ohmic and super-Ohmic case but not for sub- 
Ohmic coupling: Low energy modes become more and more important for smaller values of s in 
Eq. (|13D. 

The main price we have to pay is that our flow equations are more complicated than the 
corresponding renormalization group equations. Renormalization group equations are differential 
equations with respect to the ultraviolet cutoff and not with respect to energy differences as in 
our approach. We expect that our scheme can be useful for other models too. 
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Figure captions 

Fig. 1. Decoupling of modes as a function of the flow parameter I in our approach, < < £2 < ^3- 
J(Aoo,^) shows an algebraic decay [^. 

Fig. 2. Successive steps of the adiabatic renormalization sceme. The shaded area indicates the 
modes that are being integrated out in successive steps. 
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